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^-^ ' Excitation of multicomponent dark matter in the galactic center has been proposed as the 

source of low-energy positrons that produce the excess 511 keV 7 rays that have been ob- 

l/^ ' served by INTEGRAL. Such models have also been promoted to explain excess high-energy 

c7^ ' e observed by the PAMELA, Fermi/LAT and H.E.S.S. experiments. We investigate whether 

C^ . one model can simultaneously fit all three anomalies, in addition to further constraints from 

inverse Compton scattering by the high-energy leptons. We find models that fit both the 511 

keV and PAMELA excesses at dark matter masses M < 400 GeV, but not the Fermi lepton 

excess. The conflict arises because a more cuspy DM halo profile is needed to match the 

p^ ' observed 511 keV signal than is compatible with inverse Compton constraints at larger DM 

*~i , masses. 

C^ , 

1 Galactic cosmic ray anomalies and DM collisions 

There are several hints of unexplained sources of electrons and positrons in our galaxy, which 
could be due to collisions of dark matter (DM). The longest-standing one is the excess of 511 
keV 7 rays from the galactic center, first seen by balloon-borne detectors in the 1970's, and most 
recently measured by the SPI spectrometer aboard the INTEGRAL satellite (for a review, see 
ref.^. More recently a number of experiments have found evidence for e at higher energies, in 
excess of those understood to be coming from known sources. Among these, PAMELT^P reports 
an excess in the_positron fraction at energies of 10 — 100 GeV, while the Fermi Large Area 
Telescope (LATJk' and H.E.S.SP observe an excess of e~*" -|- e~ in the 100 — 1000 GeV energy 
range. 

Although many different astrophysical explanations have been proposed as the source of 
the low-energy positrons that produce the 511 keV signal, there is no consensus^ Pulsars have 



for example, the argument of ref.'^that low-mass x-ray binaries are most likely source has been criticized in 



been proposed as a likely source of the PAMELA and Fermi leptons (see for example refs.^^^^ 
but the uncertainties in the parameters characterizing such sources still leave room for other 
interpretations. 
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Although DM annihilations had previously been suggested as the % 
source of some of these anomalies, ref.^ was the first to point out a 
class of DM models that could potentially explain all of them (and 
a few others: the WMAP haze and the DAMA/LIBRA annual mod- 
ulation). Namely, these are models where the DM has a mass M 
near the TeV scale, and has several components that acquire natu- 
rally small mass splittings 5M ^ 1 MeV from radiative corrections. 
A new hidden sector Higgs or gauge boson with mass /^ ^ 1 GeV 
mediates annihilations of the DM into e^ but not antiprotons (since 
/i < 2mp), of which no excess has been observed by PAMELA. All of 
this can be economically achieved by assuming the hidden sector gauge symmetry is nonabelian 
and spontaneously breaks near the GeV scale. Then the mediator is one of the gauge bosons 
Bf^, which can mix with the standard model hypercharge 1^ through the dimension-5 gauge 
kinetic mixing operator A^^A'^B'^^^Y'^'^, where A" is a hidden sector Higgs field in the adjoint 
representation that gets a V EV. So me of the simplest examples involving SU(2) gauge symmetry 
were considered by us in ref.^ ^ I 



Figure 1: Inverted mass 

hierarchy for excited DM 

states. 



2 Exciting Dark Matter in the Galactic Center 

The excited DM mechanism (XDM) for explaining the 511 keV excess was first proposed in 
ref. ^^. The ground state DM particles undergo inelastic scattering to the excited state by 
XiXi ~^ X2X2; followed by decays X2 — > Xi^'^'e" into nonrelativistic e^. However a quantitative 
computation of the excitation cross section was not used there, and ref. argued that the rate 
of excitation was too small to account for the observations unless many partial waves were at 
their maxi mum values allowed by unitarity. 

In refs.l ^ ^ we numerically computed the excitation cross section by solving the Schrodinger 
equation, and showed that indeed the suspicion of ref. ^^ was correct, the rate of e^ production 
is too small, even varying all the model parameters and DM halo properties over a wide range. 
At the same time, we proposed a solution, involving the existence of a stable excited state that 
undergoes scattering X2X2 — > XsXS; followed by the decay X3 ~^ Xi^~^^~- This can have a 
smaller mass gap 5M23 ~ 100 keV which is easier to excite in DM collisions than the larger one 
(5Mi3 > 2?TT,e. This "inverted mass hierarchy" is shown in figured! 



Figure [2] (left panel) shows an example of our new 
contours for the rate of positron production compared 
to the observed rate in the M-// planel^ using the 
mass splitting SM23 = 100 keV, and the gauge coupling 
ag = 0.031 (M/TeV) required for getting the right relic 
density'MJ The DM density profile is taken to be of the 
Einasto form, 



P = Pq exp 



a 






(1) 



with Pq = 0.4 GeV/cm^ r© = 8.3 kpc, a = 0.17, 
Vg = 15 kpc, consistent with best-fit values of A^-body 

7^. 
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Figure 3: Observed angular distribution of 
INTEGRAL 511 keV signal, and theoreti- 
cal predictions for different Einasto parameter 
values Q = 0.17,0.20,0.23, r^ = 20,25 kpc. 
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Figure 2: Left: contours of the rate of positron production, log{R^+ / Roha) (for INTEGRAL 511 keV 7 rays) in 
plane of gauge boson mass fj. versus DM mass M for SM23 = 100 keV mass splitting and halo Einasto profile 
parameters a — 0.17, r^ — 15 kpc, pq — 0.4 GeV/cm'', vo — 250 km/s.^® Heavy contours match the observed 
rate. Dashed curves are contours of inverse Compton (IC) signal over IC bound. Shaded regions are excluded by 



IC constraint. Right: same, but with 5M23 = 25 keV, a = 0.20, r^ = 15 kpc, p© = 0.3 GeV/cm 



vo 



220 km/s. 
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simulations, and a high value of the circular velocity 

vq = 250 km/s. In this example, the heavy contours 

show that there exist parameters leading to a large enough rate, but these tend to disappear 

rapidly if one increases the values of 5M23 (since the X2 states do not have enough kinetic energy 

to produce X3), or a or r^ (since then p becomes too small in the central region of the galaxy, 

reducing the rate). This can be compensated by decreasing 6M on the other hand, as illustrated 

in the right panel of fig. [2l The shaded regions are ruled out by constraints on inverse Compton 

gamma rays^*^ as we discuss in the next section. Fig. [3] shows that the more cuspy DM profile 

with a = 0.17 gives a better fit to the angular distribution of the 511 keV signal. 



3 High energy e from annihilations 



Although in refs.l ^ I we showed that the XDM mechanism with inverted mass hierarchy can 
work for the 511 keV signal, we did not consider whether it could also be compatible with 
the PAMELA and Fermi lepton excesses. The same model can also explain the high energy 
leptons through annihilation to hidden sector gauge bosons, XiXi ~^ BB, followed by the 
decays B -^ e~^e~^^. However, this scenario has come under increasing pressure from various 
constraints, the most stringent being due to inverse Compton scattering of e^ on starlight in 
the galaxy, which should produce 7 rays with energies up to several hundred GeV. Demanding 
that this new source not exceed recent observations excludes the annihilating DM interpretation 
of Fermi leptons unless the galactic DM density profile is less cuspy near the centeF^ than is 
generally expected on the basis of A/^-body simulations of halo evolutionM^ This limit requires 
taking small values of 6M <, 100 keV in order to get a large enough rate for 511 keV 7 rays, as 
illustrated in fig. [2l 

The ability of the models to explain the high-energy lepton observations while respecting 
the IC constraints are summarized in figure H] taken from ref.^L^. The left figure is an example 
using a cuspy halo profile compatible with the PAMELA and 511 keV excesses, at M < 400 
GeV, while the right one shows the result of a noncuspy profile where the PAMELA and Fermi 
excesses can be marginally explained, but not the 511 keV. 
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Figure 4: Allowed regions for PAMELA and Fermi lepton excess in av-M plane, ^^ and predictions of multistate 

DM annihilation that are compatible with inverse Compton constraint. Left: for Einasto parameters a — 0.17, 

Ts = 25 kpc, pq = 0.28 GeV/cm^; Right: for a = 0.20, r^ = 30 kpc, pe = 0.3 GeV/cm^. 1// is fraction of total 

DM mass density occupied by annihilating DM ground state xi- 

4 Conclusions 

We have found that annihilating multistate DM can explain two out of three galactic cosmic ray 
anomalies, either PAMELA/Fermi or PAMELA/INTEGRAL, but not all three simultaneously. 
Although it is possible to marginally predict all the correct rates using Einasto profile parameter 
a = 0.20, the angular distribution of 511 keV 7 rays is too wide in this case. Of the two 
possibilities, the PAMELA/INTEGRAL combination seems preferable from the standpoint of 
the required DM halo parameters, since in this case we are able to adopt standard values that 
are quite compatible with A'^-body simulations of galactic structure evolution. Moreover we 
can match the anomalous lepton rates weh for PAMELA/INTEGRAL. The PAMELA/Fermi 
possibility requires stretching the halo parameters to their maximal values, while only marginally 
giving a large enough rate of leptons, yet a small enough rate of associated inverse Compton 7 
rays. 
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